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KEYWORDS Abstract In this paper, an irreversible quantum Brayton heat engine cycle model composed of two 

Finite time isomagnetic processes and two irreversible adiabatic processes is established. The quantum heat engine 

thermodynamics; works with non-interacting spin-1/2 systems and has multi-irreversibilities. By using detailed numerical 

Spin-1/2 systems; examples, this paper gives the fundamental optimal relationship of the quantum Brayton heat engine, 

Quantum Brayton heat in a general case and at high temperature limits, and analyzes the effects of internal friction and bypass 

engine cycle; heat leakage on the fundamental optimal relationship. Comparisons of the power output versus efficiency 

Power output; characteristics among this spin quantum heat Brayton engine, a quantum spin Carnot heat engine and a 

Efficiency. Brayton heat engine working with a classical working medium, are made. Three special cases, such as the 

endoreversible case, the frictionless case and the case without bypass heat leakage, are discussed. 

© 2012 Sharif University ofTechnology. Production and hosting by Elsevier B.V. All rights reserved. 


1. Introduction 

Performance analyses of heat engines have been a major 
source of thermodynamic insight, and have led to a connection 
between abstract thermodynamic theories and realizable phys¬ 
ical phenomena. Finite Time Thermodynamics (FTT), which be¬ 
gan from the study of the Newtonian law endoreversible Carnot 
engine system [1], have made considerable progress [2-8]. 
However, for some special fields and systems, such as the quan¬ 
tum amplifier, laser cooling, aerospace, superconductivity ap¬ 
plications, infra-red techniques and so on, the working medium 
of these systems obeys the law of quantum statistical mechan¬ 
ics, and the quantum characteristic of the working medium 
has to be considered in the performance investigation of these 
systems. In recent years, many authors have extended the ob¬ 
jectives of FTT to quantum thermodynamic systems, such as 
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quantum heat engines, refrigerators and heat pumps, by com¬ 
bining the semi-group approach and quantum master equation, 
and have obtained many meaningful results. 

In the performance analysis and optimization of quantum 
thermodynamic systems, the working medium are several typ¬ 
ical quantum mechanical systems, such as non-interacting har¬ 
monic oscillators [9-12] and spin-1/2 systems [13-16], ideal 
quantum gas [17-20] and microscope particles confined to a 
potential well [21], Besides the quantum Carnot cycle, quan¬ 
tum thermodynamic cycles may be the Brayton cycle, the Er¬ 
icsson cycle [14,15,22], the Otto cycle [20], the Stirling cycle 
[23-26], and so on. Similar to thermodynamic systems with a 
classical working medium, the endoreversible cycles are the 
basic models, otherwise some authors have established irre¬ 
versible models by considering the irreversibility of internal 
irreversibility [27,28], bypass heat leakage [29] and inherent 
regeneration [12,22,30], Some authors have established general 
irreversible models [15,31,32] by considering the irreversibili¬ 
ties of heat resistance, internal irreversibility (described by an 
internal irreversible factor) and bypass heat leakage, and gen¬ 
eral irreversible heat engine models [33], considering the irre¬ 
versibilities of heat resistance, internal friction (described by an 
internal friction coefficient) and bypass heat leakage. 

A Brayton heat engine cycle with a classical working medium 
consists of two adiabatic and two constant-pressure branches, 
and a quantum Brayton heat engine cycle with spin systems 
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Nomenclature 

a 

Parameter of heat reservoir (s -1 ) 

B 

Heat reservoir 

B 

External magnetic field (T) 

C e 

Dimensionless factor that describes the magni¬ 
tude of the bypass heat leakage 

c 

Parameter of heat reservoir (s -1 ) 

Es 

Internal energy of the spin-1/2 systems (J) 

H 

Hamiltonian 

h 

Reduced Planck’s constant (J • s) 

1<B 

Boltzmann constant (J/K) 

Li, L 2 , L 3 , L 4 Lagrangian functions 

M 

Magnetic moment operator 

n c 

Population of the thermal phonons of the cold 
reservoir 

P 

Power output (W) 

Q 

Amount of heat exchange (J) 


Operator in the Hilbert space of the system and 


Hermitian conjugates 

Qe 

Rate of heat flow of bypass heat leakage (W) 

Q 

Parameter of heat reservoir 

S 

Expectation value of spin operator S z 

s+,s_ 

Spin creation and annihilation operators 

S(S x ,S y ,S z ) Spin operator 

Seq 

Asymptotic value of S 

T 

Absolute temperature (K) 

t 

Time (s) 

W 

Work (J) 

Greek symbols 

a 

Intermediate variable 

p 

Temperature, /3 = 1/(k B T ) (J -1 ) 

Y+ - Y- 

Phenomenological positive coefficients 

V 

Efficiency 

X 

Parameter of the heat reservoir 

, A. 2 , A, 

3,7-4 Lagrangian multipliers 


Friction coefficient 

Pb 

Bohr magneton (J/T) 

r 

Interaction strength operator 

T 

Time (s)/cycle period (s) 

CO 

The magnetic field, a> = 2 [i B B(t) z 

Subscripts 

B 

Heat reservoir 

c 

Cold side 

h 

Hot side 

S 

Working medium system 

SB 

Interaction between heat reservoir and working 
medium system 

max, /x 

= o, C e = 0 Maximum point for endoreversible 
case 

1,2, 3,4 

Cycle states 


consists of two adiabatic and two isomagnetic field branches. 
In the investigation of quantum Brayton cycles, Feldmann 
et al. [34] were the first to establish an endoreversible quantum 
Brayton heat engine model, working with non-interacting spin 
systems, and investigated its optimal performance. Lin and 
Chen [35] established an endoreversible quantum Brayton 


heat engine model working with non-interacting harmonic 
oscillators, and investigated its optimal performance. Feldmann 
and Kosloff [36] introduced an internal friction, and established 
an irreversible quantum spin Brayton heat engine and heat 
pump models. The internal friction describes the quantum non- 
adiabatic phenomenon caused by rapid change in the external 
magnetic field in the adiabatic process, which has an obvious 
effect on the performance of the quantum heat engine and 
heat pump. Since then, the origin and effects of quantum 
friction on the performance of quantum thermodynamic cycles 
have become interesting research subjects [33,37-41], By 
taking into account heat resistance, internal irreversibility and 
bypass heat leakage, Wu et al. [32] established a generalized 
irreversible harmonic quantum Brayton heat engine model, and 
investigated the optimal performance of the two quantum heat 
engines. 

Unlike the previous study [32], this paper will establish a 
generalized irreversible quantum Brayton heat engine model 
working with non-interacting spin systems. The irreversibilities 
of heat resistance, internal friction in the adiabatic processes, 
and bypass heat leakage are taken into account in this 
model. This paper will derive the cycle period, power output 
and efficiency by using FTT, the semi-group approach and 
a quantum master equation. By using detailed numerical 
examples and illustrations, this paper will give the fundamental 
optimal relationship of the quantum Brayton heat engine and 
analyze the effects of internal friction and bypass heat leakage 
on the fundamental optimal relationship. At high temperature 
limits, the performance parameters of the quantum Brayton 
heat engine will be simplified, and the fundamental optimal 
relationship will also be obtained using detailed numerical 
examples. Comparisons of the power output versus efficiency 
characteristics between this quantum spin heat Brayton engine, 
a Brayton heat engine working with classical working medium, 
and a quantum spin Carnot heat engine, will be made. Three 
special cases (endoreversible, frictionless and without bypass 
heat leakage) will be discussed in brief. The results obtained 
are general and can enrich the FTT theory for quantum 
thermodynamic cycles. 

2. Quantum dynamics of a spin-1/2 system 

In this section, one considers a single spin-1/2 particle 
placed in a time dependent magnetic field, B, and the direction 
of the magnetic field is along the positive z axis. The 
Hamiltonian of the interaction between a magnetic moment, M, 
and the magnetic field, B, is H(t) = -M-B. For a single spin-1 /2 
system, the Hamiltonian is given by Zeng [42]: 

H s = —M ■ B = 2 |X B S • B/h = 2 \i B S z B z /h = co(t)S z /h , (1) 

where S is a spin angular momentum, which is along the 
positive M directions, p. B is the Bohr magnetron, fi is the reduced 
Planck’s constant and a>(t) = 2 p, s B(t) z . One refers to a> rather 
than B z as “the magnetic field” throughout this paper. The 
internal energy of the spin system is given by: 

E s = ^H s ) = co(s z ] j/h = ojS/h, (2) 

where: 

S = (^) = -2 tanh (?)’ (3) 
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is the expectation value of S z , p = l/(k B T), k B is the Boltzmann 
constant and T is the absolute temperature. For simplicity, one 
refers to /I rather than T as the “temperature” throughout this 
paper. 

For a spin-1/2 system coupled thermally to a heat reservoir 
(bath), the system-bath Hamiltonian is given by: 

H = H s + H sb + H b , (4) 

where H S ,H S b and H B are, respectively, Hamiltonians of 
the spin-1/2 system, the system-bath and the bath. In the 
Heisenberg picture, the motion of an operator of the spin-1/2 
system is given by the quantum master equation: 

dX if- - a 9X 

¥ = »M+i>F + 1 ° <x) - (5) 



where L D (X) is a dissipation term (the relaxation-type term 
includes the effects of H SB and H B ), which originates from a 
thermal system-bath coupling. Using a semi-group approach, 
one gives the dissipation term [43,44]: 

l d (x) = (<i + [*• Q*\ + [<i + ’ *] Qa) . ( 6 ) 

a 

where y a are phenomenological positive coefficients, and Qa 
and Q/ 1 " are system operators in Hilbert space, which are 
Hermitian conjugates. 

Substituting the system Hamiltonian, X = H s , into Eq. (5), 
gives the internal energy change rate: 



Comparing Eq. (7) with the differential form of the first law 
of thermodynamics: 

dE s dW dQ 

— = -+ —. ( 8 ) 

dt dt dt 

One can find that the instantaneous heat flow and power 


may be identified by: 

Q = (l d (Hs)) = coS/h = dQ/dt, (9) 

P = (dHs/dtj = coS/h = dW/dt. (10) 

Hence, heat and work inexact differentials are given by: 
dQ = codS/h, (11) 

dW = Sdco/h. (12) 


For a spin-1/2 system, Eq. (7) gives the time derivative form 
of the first law of thermodynamics. 

In order to obtain the change rate of spin angular momentum 
S when the spin-1/2 system couples with a heat bath, one 
chooses Q/ 1 " = S + = S x + iS y and Qa = S_ = S x — iS y in Eq. (6). 
Substituting and Qa into Eq. (5) and using the commutation 
rejections, j^S x , S y j = ihS z , |^S y , S z j = ihS x , |^S Z , S X J = ihS y and 
S 2 = S y = S 2 = ^ gives: 

S = —2h 2 (y + + y~)S — h 3 (y~ — Y+)- (13) 


Figure 1: The S — co diagram of an irreversible quantum spin Brayton heat 
engine with multi-irreversibilities. 

If the field, a>, is constant, coefficients, y + and y_, are both 
constants. Solving Eq. (13) gives: 

S(t) = S eq + [S(0) - S eq ]e _2(y++)/_)t , (14) 

where S(0) is the initial spin angular momentum, S eq = 
—h(y~ — y + )/2(y_ + y + ) is the asymptotic spin angular 
momentum, which must correspond to the equilibrium spin 
angular momentum, S eq = — |tanh(^), where is the 
temperature of the heat bath. Comparing the two expressions 
of S eq gives y-/y+ = e^“, so that one assumes simply that: 

y + = ae q ^ m , (15) 

X_ = ae^+^i 10 , (16) 

where a and q are two constants. At the weak-coupling limit, 
the explicit expressions of y + and y_ can be obtained in terms 
of correlation functions of the bath [45]. y+,y~ > 0 requires 
a > 0. If Pjw —> oo, y + —> 0 and y_ —> oo hold, it requires 
0 > q > — 1. Substituting Eqs. (15) and (16) into Eq. (13) yields 
the change rate of the spin angular momentum: 

S = -ah 2 e qf>co [2(\ + e fi,,J )S + h(e Pm - 1)]. (17) 

3. An irreversible spin quantum Brayton heat engine model 

(1) The quantum heat engine uses many non-interacting 
spin-1 /2 systems as a working medium. The spin-1/2 system in 
the engine couples mechanically to a time dependent external 
“magnetic field" and is a two level system. The direction of 
the external magnetic field is along the positive z axis. The 
magnitude of the field is not allowed to reach zero, where the 
two energy levels of the spin system are degenerate. 

(2) The quantum heat engine operates between a hot 
reservoir, B/,, and a cold reservoir, B c . The two reservoirs are 
thermal phonon systems and are at constant temperature, /S/, 
at p c , respectively. One assumes that the heat reservoirs are 
infinitely large and their internal relaxations are very strong, so 
that they can be assumed to be in thermal equilibrium. 

(3) The quantum Brayton heat engine cycle is composed of 
two isomagnetic branches (process 1 —2 with u> = m a and 
process 3^4 with m = a>b) connected by two irreversible 
adiabatic branches (processes 2^3 and 4 —1), and the S — a> 
diagram of the cycle is shown in Figure 1. 

In the two isomagnetic processes, the working medium 
couples to the heat reservoirs for a finite time, t/, and r c , 
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respectively, and exchange heat, Q 12 and Q 34 , with the heat 
reservoirs. The finite rate of the heat transfer requires p c > 
p 4 > fi\ > p 2 > Ph where P\, p 2 and p 4 are the temperatures 
of the working medium in states 1, 2 and 4, respectively. Using 
Eq. (11), one can calculate the amounts of heat the working 
medium absorbs from the hot reservoir in isomagnetic process 
1 ^ 2 : 




1 

codS — -w a (S 2 — Si) 
ft 



(18) 


From Eq. (12), one finds that there is no work done by the 
working medium along the isomagnetic processes because of 
constant a>. 

In the adiabatic processes 2^3 and 4 —> 1, there is no 
thermal coupling between the working medium and the heat 
reservoirs, so there is no heat exchange. In the adiabatic process 
2 -> 3 (or 4 -> 1), the external magnetic field varies linearly in 
time from w a (or w b ) to w b (or w a ) in a period of z a (or z b ), viz: 


w(t) = cu(O) + d>t. 


(19) 


According to the quantum adiabatic theorem [42], rapid 
change, w, of the field causes a quantum non-adiabatic 
phenomenon, which causes an increase in spin angular 
momentum, S. From Eq. (11), one understands that an increase 
in S means heat generation, i.e. work dissipation in the adiabatic 
process. One assumes that there exists a friction coefficient, /z, 
to describe the non-adiabatic phenomenon [36]: 



where t' is the time of the adiabatic process. The spin angular 
momentum is given by: 


S(t) = S(0) + ft (“) 2 t, (21) 

where 0 < t < t'.Substitutingt' = r a andt' = irrespectively, 
into Eq. (21) yields: 

S 2 = S(T a ) = S 3 - ft/x 2 / t q , (22) 

S 4 = S( T b ) = Si - h/i 2 /z b , (23) 

where S, = -§ tanh ^ and S 3 = tanh ^ are the 
spin angular momentums at states 1 and 3, respectively. From 
Eqs. (3), (22) and (23), one can obtain the temperatures of the 
working medium at states 2 and 4, respectively: 


P2 

Pa 


— tanh 1 (1 

V 


2 

U>1 j 


— tanh -1 h 


,a„h ft "* + 2 ' i2 ), 

(24) 

tanl/'"- + 2, ‘ 2 ). 

(25) 


From Eqs. (3), (24) and (25 ), one can obtain the upper bound 
of P\ and the lower bound of p 3 : 


2 / Pc<*>b 2u 2 \ 

Px < p :r = — tanh 1 (tanh —-I , (26) 

(V a \ 2 T b J 

2 _1 / PhU>a 2/Z 2 \ 

Pi > Pir = — tanh (tanh -) . (27) 

(Db \ 2 T a J 


The work done by the working medium along the adiabatic 
processes is equal to the internal energy changes, due to the fact 
that there is no heat exchange between the working medium 
and heat reservoirs. Using Eqs. (7), (19) and (21 ), one can obtain: 


Wr 


W 4 i 


P ta 1 P r a 1 P t 

= - dEs = — — / Sdw — - I 
Jo h Jo ft Jo 


C r a 1 P ta 

Sdw -/ wdS 

ft . 

/S 3 /z 2 \ LL 2 (w a + CO b ) 

(o, c - Wbl | — j - . 


r r b i f Tb 1 r r b 

~J. dCs = -nl Sd -»i ' 


wdS 


(w b - w a ) 


2r b 


H z {w a + w b ) 
2z b 


(28) 


(29) 


The second parts of W 23 and W 41 , i.e. M anc i ^ 

ZTq ZTb 

is the work done against the internal friction along the adiabatic 
processes, 2^3 and 4 —► 1, respectively. 

(4) The quantum Brayton heat engine has a bypass heat 
leakage between hot and cold reservoirs. It arises from the 
thermal coupling action between the hot and cold reservoirs by 
the working medium. 

As mentioned above, the hot and cold reservoirs are both 
thermal phonon systems, B b and B c , respectively. One assumes 
that the creation and annihilation operators of the phonons in 
the hot and cold reservoirs are, respectively, fa ( “, ft+ and b~, 
and the phonon frequency and population of the cold reservoir 
are, respectively, w c (t) and n c = \/(e hWc ^ c — 1). When there 
exists a thermal coupling action between the hot and cold 
reservoirs by the working medium, the heat reservoirs become 
an open system. Similar to the derivation of S in Section 2, one 
can derive the time derivative of the population of the cold heat 
reservoir, n c , by solving the quantum master equation of the 
heat reservoirs [29]: 


h c = -2 _ 1)nc _ 1 ] (3 0 ) 

where c and X are constants. By using the instantaneous heat 
flow definition for the harmonic oscillators system, Q = 
wh [ 10 ], the rate of bypass heat leakage (rate of heat flow from 
the hot reservoir to the cold reservoir) is given by: 


4 . = C e hw c h c = 2C e chw c e^ hmc [l - (e hfth<0c - l)n c ], (31) 

where C e is a dimensionless factor one introduces to describe 
the magnitude of the bypass heat leakage. One assumes that 4 
is a constant, therefore, the bypass heat leakage quantity, Qe, 
per cycle, is given by: 


(^ = 41 = 2C e cftw c e A,lft “ c [l - (e nph ° lc - l)n c ]r. (32) 

The model studied in this paper is similar to the model of 
classical generalized irreversible simple Brayton heat engines 
with multi-irreversibilities, such as heat resistance, internal 
irreversibility and bypass heat leakage [46-51], 


4. Cycle period 


From Eq. (17), one can calculate the time of isomagnetic 
processes 1^2 and 3^4, respectively: 


Tf. = 


C 2 ds _ 

~ is, J " 


1 


x In 


2 a ft 2 (ePb ma + 1) 

tanh(/l,,a) a / 2 ) — tanh(^io> a / 2 ) 


tanh( yS/ 1 <x> a / 2 ) — tanh(y 0 3 a>i,/ 2 ) — 2jj?/z a 


, (33) 
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L 


s4 dS 
S3 J 

x In 


1 


2 a h 2 e q P cC ° b (eP c0,b + 1 ) 
tanh( / S c a)/,/2) — tanh(/l 3 a> b /2) 


_tanh(/3 c ai(,/2) — tanh(/lia> a / 2 ) — 2 /z 2 /T h _ 
The cycle period is given by: 

r = T h + T c + T a + T b . 

5. Fundamental optimal relationship 


Using Eqs. (28) and (29) yields the total work output per 
cycle of the quantum Brayton heat engine: 




l.(K 

(34) 

II 

u 

0 

0.8. 


II 

0.6. 

(35) 

cC 

0.4. 


Sr 

0.2. 



0.0 S 

0.3 


W = (f) dW = W U + W 23 + W 34 + W 41 

' P\O) 0 


1 

-(a) a - co b ) 


tanh 


— tanh 


PiU>b 



0.3 


Figure 2: 
tures” 


0-8 0.6 


The dimensionless power output P/P m ax,Ai=o,c e =o versus “tempera- 
and /} 3 . 


-/X 


CO a + (Ob 
Ta T b 


(36) 


Using Eq. (35) with Eq. (36) yields the power output of the 
quantum heat engine: 


P = 


W 

T 



U>b) 




- M 


2 



(37) 


Using Eqs. (18), (32) and (36) yields the efficiency of the 
quantum heat engine: Eq. (38) is given in Box I, where Q;, = 
Q .12 + Qe is the total heat released by the hot reservoir per cycle. 

Using numerical calculation, one can plot three-dimensional 
diagrams of the power output (P/Pmax,/i=o,c e =o, > A) and ef- 
ficiency (??, p u p 3 ) from Eqs. (37) and (38) where P m ax,^=o.c e =o 
is the maximum power output for the endoreversible case (/z = 
0, C e = 0), as shown in Figures 2 and 3. According to Refs. 
[34,36], h = 1 and k B = 1 are set in the numerical calcula¬ 
tions for simplicity. The parameters used in the numerical cal¬ 
culations are a = c = 2, q = X = —0.5, f} b = 0.2, = 

l,r a = T b = 0.01, a)] =5 ,a >3 = 3,a> c = 0.2, fj. = 0.002 
and C e = 0.05. From Figure 2, one can see that there exist 
optimal temperatures, fi\ and f} 3 , which lead to the maximum 
power output of the quantum Brayton heat engine. Affected by 
the internal friction and bypass heat leakage, the maximum di¬ 
mensionless power output (P/P max , M =o,Q=o)max < 1 . Figure 3 
shows that there also exist optimal temperatures, P\ and f} 3 , 
which lead to maximum efficiency when there exists a bypass 
heat leakage. 

For the quantum Brayton heat engine, to determine the 
maximum power output for fixed efficiency or maximum 
efficiency for a fixed power output, one introduces Lagrangian 
functions; Li = P + and L 2 = rj + A 2 P, where Xj and 
X 2 are two Lagrangian multipliers. Theoretically, by combining 
Eqs. (37) and (38) and the extreme conditions, one can obtain 
the optimal relation between fii and /} 3 , and then obtain the 
fundamental optimal relationship. The extreme conditions are 
given by the Euler-Lagrange equation: 

dLi/dPi = 0, dLi/dp 3 = 0, (39) 


0.400, 


0.395, 


R- 0.390, 



0.3 0.3 


Figure 3: The efficiency r] versus “temperatures” pi and 


However, because of the strong complexity and nonlinearity 
of these equations, one cannot solve these equations analyt¬ 
ically, and obtain fundamental optimal relations analytically. 
With the same parameters used in the calculation of Figure 2, 
one can plot the dimensionless power output versus efficiency 
characteristic curves, P/P mlx .ii=o.c e =o — r h by solving the Eu¬ 
ler-Lagrange 

Eq. (39) or (40), numerically, as shown in Figure 4. One can see 
from Figure 4 that the efficiency is a constant (for given co b /oj a ) 
when there exists no internal friction and bypass heat leak¬ 
age. When a bypass heat leakage exists, the P/Pmax,/*=o,Q=o — V 
curve is monotonic increasing. When internal friction exists, the 
P/Pmax,ii=o,Ce=o — d curve is parabolic-like. When both internal 
friction and bypass heat leakage exists, the P/P max ,j/=o,c e =o — 
ri curve becomes loop-shaped. The internal friction decreases 
both power output and efficiency, while the bypass heat leak¬ 
age only decreases efficiency. 

6. Fundamental optimal relationship at high temperature 
limit 

When the temperatures of heat reservoirs and the working 
medium are high enough, i.e. pw 1 , the results obtained 
above can be simplified. For example, Eqs. (3), (31) and 
(33)-(38) can be simplified at the first order approximation: 
Eqs. (41)-(48) are given in Box II where: 


2 cfia) c (l + XhPljCDc) 
(40) “ = J c 


3L 2 /dPi = 0 , 


3 L 2 /dp 3 = 0, 
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w 

ft = — 


K - co b ) [tanh (&j±) - tanh (^)] - 2 m 2 (f + 

Qi> co a [tanh — tanh _ iitij _|_ 4C £ ,cfl&) c e ;lfi ^' :i ' c [l — (e n ^ a,c — l)n c ]t 


(38) 


Box I: 


s = --?*, 

2chw c n + xhp b a> c ) 

Qe ~ C e - „- ~ Ph) = CMPc ~ Ph), 


Th 


4 ah 2 


In 


TaOLaiPh P l) 


TaOJaPh PiO^b T a 4// 2 

Tb<*>b(Pc ~ Pi) 

4a ft 2 \_T b co b p c ~ P\co a Tb ~ 4/x 2 


1 


In 


1 


In 


0) a (x>bT a Tb(Ph P\)( Pc Pi) 


4a ft 2 (z a w a p b - p 3 CL> b T a - 4ii 2 )(T b ObPc - P\CO a Tb - 4m 2 ) 

... (ft) a - C0 b )(P 1 0) a - Pi(O b ) 2 (, M b 

w = ---Ml- 1 - 

4 V T, 


+ T a + T b , 


Tb 

2 .,2 


P = 


ft = 


a ft (a> a — w b ){fi\u> a — Pi(D b ) — 4aft fi ( T a a> b + T b w a ) 

T a T b In - - - “ aa>b f ° tbi K a C ~ a 3 * --pTr + 4<3 ft 2 T a X b (z a + T b ) 

(t a iu„- ft r 0 - 4/i 2 Kq mj, ^, m 0 r 6 - 4/j 2 ) 

a ft 2 z a z b (co a - u>b){P\C 0 a - p 3 co b ) - 4a ft 2 /z 2 (z a co b + z b co a ) 


a ft 2 o) a z b {z a u> a Pi - p 3 z a co b - 4m 2 ) + z a z b C e a(p c - Ph) 


In 


u> a u> b z a z b (p b Pl)(Pc Pi) 


(z a w a p b - p 3 co b x a - 4ix 2 )(x b a>bPc ~ P\u> a z b - 4m 2 ) 


+ 4a fti 2 (x a + z b ) 


(41) 

(42) 

(43) 

(44) 

(45) 

(46) 

(47) 

(48) 


Box II: 



Figure 4: The effects of n and C e on dimensionless power output P/ 
Pm&,n=o,C'=o versus efficiency 

Similar to Figures 2 and 3, one can plot three-dimensional 
diagrams of the power output (P/P miXtfl= o,c e =o, Pi, Pi) and ef¬ 
ficiency (r), Pi, p 3 ) from Eqs. (47) and (48), where P max ,^=o,c £ =o 
is the maximum power output for an endoreversible case (m = 
0, C e = 0) at high temperature limits. For simplicity, these two, 
three-dimensional diagrams are not given here. At high temper¬ 
ature limits, one can see that the relationship between dimen¬ 
sionless power output (or efficiency) and temperatures, Pi and 
p 3 , is similar to that of a general case. The dimensionless power 


output has a maximum, and when bypass heat leakage exists, 
the efficiency has a maximum that has a nonzero correspond¬ 
ing power output. 

For the quantum Brayton heat engine, to determine the 
maximum power output for fixed or maximum efficiency for a 
fixed power output, one introduces Lagrangian functions, L 3 = 
P + X 3 q and L 4 = q + X 4 P where and X 4 are two Lagrangian 
multipliers. The Euler-Lagrange equations become: 


dLi/dPi 

= 0, 

dL 3 /dp 3 = 0, 

(49) 

or: 

dU/dPi 

= o, 

dL 4 /dp 3 = 0. 

(50) 


With the help of Eqs. (47) and (48), one can solve 
Eq. (49) or (50), numerically, and obtain the fundamental 
optimal relationship of the quantum Brayton heat engine at 
high temperature limit, as shown in Figure 5. According to 
Refs. [34,36], the parameter values used in the numerical 
calculations are: 

a = c = 2, q = X = —0.5, p b = 0.001, 

Pc = 1/290, z a = z b = 0.01, a)i = 10, 

a> 3 = 8 , a> c = 6 . 

From Figure 5, one can see that the efficiency, q, is constant 
(for given o) b /a) a ) and the P/P ma x, AI =o,c c =o - ft curve is a 
beeline when no internal friction and bypass heat leakage 
exist. The internal friction makes the P/P max #(=0 ,c E =o — ft curve 
become parabolic-like, while the bypass heat leakage makes 
the P/P max , # i=o,c e =o — ft curve become monotonic increasing. 
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Figure 5: The Effects of fi and C e on dimensionless power output P/ 
P ma x,M=o,c e =o versus efficiency r/ at high temperature limit. 

When both internal friction and bypass heat leakage exist, the 
P/Pnax,n=o,Ce=o ~ P curve becomes loop-shaped. 

The fundamental optimal relationship of the irreversible 
spin quantum Brayton heat engine is different to those of 
an irreversible spin quantum Carnot heat engine [31] and 
an irreversible Brayton heat engine with classical working 
medium [51], The P/P m ax, / u=o,c e =o — b curves of the spin 
quantum Carnot heat engine and the Brayton heat engine 
with classical working medium are both parabolic-like for 
an endoreversible case, and the bypass heat leakage changes 
the P/P m ax,/i=o,Q=o — b curves to loop-shaped. The internal 
irreversibility only affects the power output versus efficiency 
characteristic, quantitatively. The difference means that the 
quantum characteristic of the working medium does affect the 
performance of the quantum heat engine, and it is necessary 
for one to consider the quantum characteristic of the working 
medium in the performance investigation of a quantum heat 
engine. 


At high temperature limits, the cycle period and power 
output can be simplified to: 


1 COaCObtfk - ft)(ft - ft) 

r = -=■ In-, 

4aft Oaft - fta>/,)(a>i,ft - fta>„) 

P = a h 2 (a> a — U>b){P\U>a ~ ft®/)) 

co a a>b(fih ~ P\)(Pc ~ ft) 

x In- 

. ((Waft - P 3 0)b)(o>bPc - fta>a)_ 


(54) 


(55) 


Case 2: Frictionless case (i.e. /x = 0 and C e ± 0). The 
irreversibilities in the cycle are heat resistance and bypass heat 
leakage. Similar to the endoreversible case, by neglecting the 
time of the adiabatic processes (i.e. r a = rb = 0), one can 
simplify the efficiency as: Eq. (56) is given in Box III. 

The cycle period and power output are independent of the 
bypass heat leakage, therefore, the cycle period and power 
output are still Eqs. (51) and (52), respectively. Based on 
Eqs. (52) and (56), one cannot derive the fundamental optimal 
relation analytically. Using numerical calculation, one obtains 
the P/P m ax,/j=o,Q=o — b characteristic curve, as shown in 
Figure 4 (line 3), and the curve is monotonic increasing. 

At high temperature limits, one can simplify Eq. (56) as: 
Eq. (57) is given in Box IV. 

From Eqs. (55) and (57), one can derive the fundamental 
optimal relation, analytically, of the quantum Brayton heat 
engine in a frictionless case as: 


P = C e a(ft — ft)(<w a — co b )——-ft-. (58) 

M a (l — rf) — a>b 

Figure 5 (line 3) gives the P/P max ^=o,c e =o ~ V curve in the 
frictionless case and the curve is monotonic increasing. 

Case 3: The case without bypass heat leakage (i.e. /x ^ 0 and 
C e = 0). The irreversibilities in the cycle are heat resistance and 
internal friction. In this case, one can simplify the efficiency as: 


*1 = 



M a \ 


1 


7. Three special cases 


The results obtained above include the fundamental optimal 
relations in an endoreversible case, a frictionless case and a case 
without bypass heat leakage. 

Case 1: Endoreversible case (i.e. /x = 0 and C e = 0). The 
solo irreversibility in the cycle is heat resistance. The time of 
the two adiabatic processes is negligible (i.e. r a = z h = 0) 
compared to the time spent on the isomagnetic processes. The 
cycle period, power output and efficiency can be, respectively, 
expressed simply as: 


1 

2a h 2 e^ hU ‘ a (e^ h01a + 1) 

tanh(fto> a /2) - tanh(ft<w fl /2) 

tanh(fto> a /2) — tanh(fta>/,/2) 

1 

H-~- 

2 a fr e q P cC ° b (e^ cC ° b + 1) 

1 "tanh(fto>b/2) — tanh(fto)f,/2) 
_tanh(fta> b /2) — tanh(fta> a /2) 


(51) 


XX _ \ U U 1 

tanh(fto> a /2) — tanh(ft&>/,/2) — 2/x 2 /r a 

The cycle period and power output are independent of 
bypass heat leakage, therefore, the cycle period and power 
output are still Eqs. (35) and (37), respectively. Based on 
Eqs. (37) and (59), one cannot derive the fundamental optimal 
relation analytically. Using numerical calculation, one obtains 
theP/P max>M=0 ,c e =o — h characteristic curve of the quantum heat 
engine without bypass heat leakage, as shown in Figure 4 (line 
2), and the curve is parabolic-like. 

At high temperature limits, one can simplify Eq. (59) as: 


\ 4/ft 

---5-- (60) 

) (T7i« a ft - r a w b p 3 - 4/x 2 ) 

Based on Eqs. (47) and (60), one cannot obtain the 
fundamental optimal relation of the quantum Brayton heat 
engine analytically. Using numerical calculations, one obtains 
the P / F m ax.,j=o.c e =o — rj characteristic curve, as shown in 
Figure 5 (line 2), and the P/P max ,/j.=o,q=o ~ *1 curve is parabolic¬ 
like. 


co b u> b [ r a 

i) = \ -( 1-1 - 

COq V T b 


p 


1 

2 (ft, ° 


Mb) 




(52) 


8. Conclusions 


. M b 

]]= l -. 

M a 


(53) 


For the irreversible quantum spin Brayton heat engine, 
the numerical examples show that both power output and 
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(w a - ft»(,)[tanh( / Si«„/2) - tanh(y0 3 aJt,/2)] 

&> Q [tanh(/Si&> a /2) — tanh(j6 3 a>f,/2)] + 4C e chco c e xh ^ mc [l — — l)n c ]r ' 

(56) 

Box III: 

a h 2 (u> a — o)b)(fi\ 0 ) a — ftcot) 
a ft 2 w a (P\(D a M + CMPc A) In 

(57) 

Box IV: 


efficiency have maxima. The power output versus efficiency 
characteristic curve is a beeline, and the efficiency is constant 
in the endoreversible case. Internal friction makes the power 
output versus efficiency characteristic curve a parabolic-like 
line, while the bypass heat leakage makes it a monotonic 
increasing line. When both internal friction and bypass heat 
leakage exists, the curve becomes a loop-shaped line. At 
high temperature limit, the power output versus efficiency 
characteristics is similar to that in a general case, and there 
is only a quantitative difference between them. Comparisons 
of the power output versus efficiency characteristics between 
this spin quantum heat Brayton engine, a Brayton heat engine 
working with classical working medium and a spin quantum 
Carnot heat engine, are made. 

The obtained results are general; they can enrich the FTT 
theory for quantum thermodynamic cycles and can provide 
some guidelines for the optimal design of a real quantum heat 
engine. 
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